Linear- Temperature Dependence of Static Magnetic Susceptibility in LaFeAsO from 

Dynamical Mean-Field Theory 
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In this Letter we report the LDA+DMFT (method combining Local Density Approximation with 
Dynamical Mean-Field Theory) results for magnetic properties of parent superconductor LaFeAsO 
in paramagnetic phase. Calculated uniform magnetic susceptibility shows linear dependence at 
intermediate temperatures in agreement with experimental data. For high temperatures (>1000 
K) calculations show saturation and then susceptibility decreases with temperature. Contribu- 
tions to temperature dependence of the uniform susceptibility are strongly orbitally dependent. 
It is related to the form of the orbitally-resolved spectral functions near the Fermi energy with 
strong temperature dependent narrow peaks for some of the orbitals. Our results demonstrate that 
linear-temperature dependence of static magnetic susceptibility in pnictide superconductors can be 
reproduced without invoking antiferromagnetic fluctuations. 

PACS numbers: 71.27.+a, 71.10.-w, 74.70.Xa 
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Introduction. LaFeAsO is the first and the most 
investigated representative of novel iron pnictide su- 
perconductors. All known to date compounds of this 
class: ReFeAs0i_2;F:r {with Re=La,Ce,Nd,Sm,Gd) and 
AFe2As2 (A=Ca,Sr,Ba) contain FeAs layers and exhibit 
a common phase diagram with antiferromagnetic spin- 
density wave (SDW) order appearing below Tn for parent 
compounds. The antiferromagnetism is suppressed by 
doping or pressure with simultaneous appearance of su- 
perconductivity. That picture resembles cuprates where 
superconductivity is associated with suppression of mag- 
netism by electron or hole doping, although the undoped 
parent compounds for the FeAs-superconductors are not 
Mott-Hubbard insulators but multiband metals. These 
facts suggest the complex interplay of magnetism and 
superconductivity in iron pnictide materials. 

Recently it was found that magnetic properties of these 
materials show anomalous behavior even in paramagnetic 
phase. A universal linear-temperature dependence of the 
static magnetic susceptibility was observed [2!-^ in both 
parent and superconducting compounds in paramagnetic 
normal state above antiferromagnetic Tn or supercon- 
ducting Tc transition temperatures. This non-Pauli and 
non-Curie- Weiss-like dependence cannot be understood 
within standard theories of magnetism. The authors of 
@i propose that this anomalous behavior appears due 
to the presence of a short-range antiferromagnetic fluctu- 
ations in the paramagnetic phase of iron pnictides. The 
antifferomagnetic fluctuations were also argued to en- 
hance the nonanalytic (linear) term in the temperature 
dependence of susceptibility in a two-dimensional Fermi 
liquid 0, giving a possibility to provide a contribution 
comparable with experimental observations. 

An important issue is the effect of the Coulomb in- 
teraction on the electronic and magnetic properties of 
iron pnictides It is generally accepted now that 



Coulomb correlation effects in iron pnictides are not as 
strong as in cuprates where parent compounds are wide 
gap Mott insulators. Experimental spectra do not show 
features corresponding to Hubbard bands that are typ- 
ical for strongly correlated materials. However moder- 
ate correlations lead to a significant renormalization of 
the electronic states near the Fermi energy with effec- 
tive electronic mass value m* jm w 2 and corresponding 
band narrowing of the LDA band structure is observed 
in ARPES 

In this Letter we report the first ab-initio computa- 
tional results for uniform magnetic susceptibility of par- 
ent superconductor LaFeAsO in paramagnetic phase ob- 
tained with LDA+DMFT, which takes into account the 
effect of the on-site correlations. A linear increase of mag- 
netic susceptibility is found below 1000 K with satura- 
tion and decreasing for higher temperatures. As DMFT 
takes into account only local correlations, the anomaly in 
temperature dependence of susceptibility observed in our 
calculations is not connected with antiferromagnetic fluc- 
tuations. It is shown that the observed temperature de- 
pendence of susceptibility can be associated with narrow 
temperature dependent peaks in orbitally-resolved densi- 
ties of states, arising due to local correlations and located 
«100 meV below the Fermi energy. These peaks yield in- 
crease of the susceptibility with temperature, similar to 
that in metamagnetic systems. In the high-temperature 
regime the electrons loose their coherence and the Curie- 
Weiss behavior of susceptibility is restored. 

Method- The LDA-I-DMFT scheme proceeds in two 
steps: (i) construction of the effective Hamiltonian from 
converged LDA calculation and (ii) solution of the corre- 
sponding DMFT equations. In the present work the pro- 
jection procedure onto Wannier functions 12| was used to 
obtain an effective 22-band Hamiltonian i?LDA(k) which 
incorporates five Fe d, three O p and three As p orbitals 
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per formula unit (two formula units per unit cell). 

The DMFT self-consistent equations were solved 
for imaginary Matsubara frequencies. The effective 
impurity problem was solved by hybridization func- 
tion expansion continuous-time quantum Monte-Carlo 
method (CT-QMC)[ll with Coulomb interaction taken 
in density-density form. The interaction matrix U^^, 
was parametrized by parameters U and J according to 
procedure described in . In the present work we used 
[7=4 eV and J=l eV obtained via constrained LDA cal- 
culations procedure in Ref. [l^. Calculations were per- 
formed in the paramagnetic state at the inverse temper- 
atures 13 =4^30 eV-i (from 387 to 2900 K). The en- 
ergy dependence of the self-energy T,{uj) on the real axis 
needed to calculate spectral functions was obtained by 
the analytical continuation via Fade approximants [lal 
according to procedure described in Ref. [l^]. 

Orbitally resolved DMFT densities of states were cal- 
culated as 

Mio) = --Im^[(c. + - H{k) + ±{uj)]T^, (1) 



where ^f(k) is the effective Hamiltonian and fi is the self- 
consistent chemical potential. The Hamiltonian -ff (k) is 
obtained by subtracting a Coulomb interaction energy 
Edc already present in LDA (so-called double-counting 
correction) from d-d block of ffLDA(k). The double- 
counting has the form Edc = [^("-dmft — 5) where 
nriMFT is the total number of 3rf electrons obtained 
within DMFT in the absence of the external field (such 
approximation is justified for small fields used in uniform 
susceptibility calculation) and U is the average Coulomb 
interaction parameter for the d shell. This choice of a 
double-counting term is not unique and other variants 
can also be used [3] ■ The form of Edc used in the present 
work yields good results for transition oxide compounds 

13. 



[19| and pnictide materials [15|, [17 



The uniform magnetic susceptibility was calculated as 
dMiT) 5K(r) - n^iT)] 



X{T) 



dEh 



dEh 



(2) 



where M{T) is the field-induced magnetization and E^ is 
the energy correction corresponding to the applied field. 
For all considered temperatures the absence of polariza- 
tion at zero field was checked. The derivative was com- 
puted as a ratio of the magnetization M(T) and the en- 
ergy Eh at small fields (where M{T) is a linear function 
of Eh). To minimize numerical error due to Monte- Carlo 
procedure for each T a series of x{T) calculations corre- 
sponding to different magnetic fields was done. 

Results and discussion - In FiglT] we show calculated 
x(T) in comparison with the experimental data of Klin- 
geler et al.^. For temperatures below wlOOO K the cal- 
culated x{T) demonstrates roughly linear increase. To 
emphasize linear character of the susceptibility curve in 



low-temperature region the least-square fit by a straight 
line to the obtained data is also shown in Fig[TJ For 
T > 1000 K the susceptibility shows saturation and then 
decreases with temperature approaching to Curie- Weiss- 
like behavior at very high temperatures. To clarify the 
origin of the linear dependence of x(2^) iu the low tem- 
perature regime it is useful to calculate orbitally resolved 
susceptibilities, Xm = dMm/dEh- Different d-orbitals of 
iron give diverse contributions to calculated x{T) (see up- 
per panel of Figl2]). The contribution of 3z^ — orbital 
posesses pronounced largest positive slope and provides 
the major source for linear increase of the total suscepti- 
bility below 1000 K. 

It can be shown that the observed temperature de- 
pendence of XmiT) originates from pecularities of or- 
bitally resolved densities of states Ai{uj), obtained within 
LDA-t-DMFT (see Fig. E]). The calculated functions 
Ai{uj) are in good agreement with previously published 
results of Aichhorn et al. 2l| . The correlations result in 
renormalization of the LDA bands with corresponding 
enhancement of the quasiparticle mass m*/m k,2 and 
lead to appearance of peaks near the Fermi energy. The 
peak corresponding to iz^ — orbital is especially sharp 
and narrow while the peaks for other orbitals are much 
more broad. The peaks «100 meV below the Fermi en- 
ergy are absent in the LDA band structure and originate 
purely from the local correlations; they display strong 
temperature dependence, increasing in amplitude, nar- 
rowing and approaching the Fermi level with decreasing 
T. 

These pecularities of orbitally resolved densities of 
states allow to explain qualitatively the anomalous be- 
havior of the susceptibility. The origin of increase of 
x(T) with temperature is similar to that in systems with 
van-Hove singularities near the Fermi level 22 - ij] , where 
increasing temperature yields activation of the electronic 
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FIG. 1: (Color online) Uniform susceptibility x(2^) of 
LaFeAsO calculated within LDA-I-DMFT (squares) in com- 
parison with experimental data of Klingeler et aZ.[3| (circles). 
The straight line corresponds to least-square fit of the region 
below 1000 K. The inset shows the full experimental suscep- 
tibility curve. 
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states with the energy of the peak position, which results 
in increase of magnetic susceptibiUty. When the energy 
ksT {T > 1000 K) becomes larger than distance from the 
peak to the Fermi level (100 meV), the linear increase of 
uniform susceptibility stops and the curve goes to sat- 
uration. The same mechanism applies also to the con- 
tribution of — 2/^ orbital, which does not show peak 
of the spectral function, but a minimum very near the 
Fermi level and increase at larger energies. 

To verify that the observed pecularities of the single- 
particle properties are responsible for non-monotonic 
temperature dependence of x{T), we have estimated the 
uniform susceptibility as a convolution of the DMFT in- 
teracting Green's functions neglecting vertex corrections. 



X 



(3) 



where G,^„'(k,icj) = [{uj + n)i - H{k) + ±{uj)],J^^,. 
Temperature dependence of the calculated susceptibili- 
ties Xm{T) is shown in the lower panel of Figl2] The 
obtained curves qualitatively reproduce all features of 
the directly calculated susceptibilities of — and 
3z^ — l orbitals, including their slope at low temperatures 
(T<1000 K). Note that the off-diagonal components of 
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FIG. 2: (Color online) Orbitally resolved Fe 3rf susceptibilities 
Xm vs temperature of LaFeAsO from the LDA-I-DMFT cal- 
culation obtained as derivatives of the magnetization(upper 
panel) and estimated via convolution of the corresponding 
Green's functions (lower panel). 



TABLE I: Imaginary part of the self-energy (eV) at the Fermi 
energy Ep at different temperatures in LaFeAsO. 



Temperature 


387 K 580 K 


1160 K 


ImE,j;(Ep) 


-0.142 


-0.242 


-0.454 


ImSy2,a;z(Ep) 


-0.131 


-0.163 


-0.306 


ImE322_^2 (Ep) 


-0.054 


-0.092 


-0.228 


ImE3.2 _^2 (Ep) 


-0.053 


-0.101 


-0.334 



the Green functions also contribute partly to the linear 
temperature dependence of susceptibilities. For xy^ xz 
and yz orbitals the vertex corrections appear to be more 
important. 

The importance of the vertex corrections for some or- 
bitals can be related to the loss of the electron coherence 
caused by electronic correlations, which can be traced 
from the electronic self-energy (Table 1). For example 
the imaginary part Ti^z.yz at the Fermi level is 131 meV 
at the lowest considered temperature and increases to 
306 meV at T=1160 K. The loss of the electron coher- 
ence is much less pronounced for 3z^ — and x^ — y^ 
orbitals where ImS(EF) are three times smaller than the 
corresponding values for other orbitals at the lowest tem- 
perature. 

Due to loss of coherence of electronic excitations in 
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FIG. 3: (Color online) Temperature evohition of LaFeAsO 
Fe 3rf spectral functions in vicinity of the Fermi energy (0 
eV) calculated within LDA-I-DMFT. 
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part of the bands, the local susceptibility behaves very 
differently from the uniform susceptibility. In Fig. |4] 
we plot temperature dependence of the Fe-d local spin 
susceptibility XiodT) and orbital contributions to its in- 
verse value. The calculated Xioc{T) decreases monoton- 
ically with temperature in agreement with the results 
reported in . We have also calculated the mean square 
of the local moment < ml > as a function of tempera- 
ture, which is almost temperature independent increasing 
approximately by only 10% of its value in the tempera- 
ture region T=387 1160 K. Analysis of the orbitally- 
resolved inverse susceptibilities (see inset of Fig. |4|) shows 
that the contribution of xz and yz orbitals is almost lin- 
ear in temperature, which goes along with loss of elec- 
tron coherence of the corresponding states, providing a 
possibility of the local moment formation in these or- 
bitals. At the same time, the electronic states origi- 
nating from 3z^ — orbital remain more quasiparti- 
cle {lm.'Exz,yz(EF) ~Slm.'E3z^_r2(Ep)). Such a picture 
is reminiscent of the proximity to the orbital-selective 
Mott transition, proposed recently for iron [25|, iron 



pnictides 26| and iron oxide[27]. Quasiparticle nature 
of 3z^ — states allows to justify the above discussed 
scenario of the linear dependence of susceptibility due 
to contribution of peculiarities of the density of states 
in — band (narrow temperature sensitive peak in 
spectral functions at 100 mcV below the Fermi level). 

Conclusion- By employing the LDA-HDMFT(CT- 
QMC) method we computed paramagnetic uniform 
magnetic susceptibility of the parent superconductor 
LaFeAsO. The obtained curve reproduces linear behav- 
ior at low temperatures observed in experimental data of 
Klingeler et al.^. We argue that the low-temperature 
linear increase of susceptibility in LaFeAsO comes from 
the presence of a sharp temperature dependent peak in 
the spectral function at 100 meV below the Fermi energy. 
Our results demonstrate that increase of the susceptibil- 
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FIG. 4: The total local spin susceptibility versus tempera- 
ture and temperature dependence of orbitally resolved in- 
verse of the local spin susceptibility(inset) obtained within 
LDA-^DMFT. 



ity can be understood within single-site dynamical mean- 
field approach neglecting spatial magnetic fluctuations. 
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